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Abstract — This paper is concerned with optimum diversity re- 
ceiver structure and its performance analysis of differential phase 
shift keying (DPSK) with differential detection over nonselective, 
independent, nonidentically distributed, Rayleigh fading chan- 
nels. The fading process in each branch is assumed to have an 
arbitrary Doppler spectrum with arbitrary Doppler bandwidth, 
but to have distinct, asymmetric fading power spectral density 
characteristic. Using 8-DPSK as an example, the average bit error 
probability (BEP) of the optimum diversity receiver is obtained 
by calculating the BEP for each of the three individual bits. 
The BEP results derived are given in exact, explicit, closed-form 
expressions which show clearly the behavior of the performance 
as a function of various system parameters. 

I. INTRODUCTION 

The receiver structure and bit error probability (BEP) per- 
formance of differential phase shift keying (DPSK) with differ- 
ential detection over nonselective, independent and identically 
distributed (i.i.d.), Rayleigh fading channels with combining 
diversity reception have been well known in the literature 
[1]— [4]. However, reaserch shows that in some practical 
systems, the independent, non-identically distributed (i.n.i.d.) 
channel model is more accurate [5], [6]. In i.n.i.d. channel, 
the fading processes and possibly the additive, white Gaussian 
noise (AWGN) on the diversity branches have non-uniform 
power profiles which are distinct from one another. The effect 
of the nonidentical diversity branch statistics on the receiver 
structure is studied in [7]. Recently, based on the maximum a 
posteriori probability (MAP) criterion, an explicit structure of 
the optimum combining differential receiver and a complete 
set of closed-form BEP expressions and their Chernoff upper 
bounds, for 2-, 4- and 8-DPSK, both with optimum combining 
reception and suboptimum combining reception, are derived 
in [8] — [10]. The purpose of this paper is to provide a further 
extension. The results derived in this paper, together with those 
in [8] — [10], form a benchmark counterpart to the classic ones 
for the i.i.d. channel given in [1] — [4]. 

In a Rayleigh channel, the fading gain is usually modeled 
as a zero-mean, stationary, complex, Gaussian random pro- 
cess. The most widely accepted model [1] — [10] is that the 
spectrum of the fading process over each diversity branch is 
symmetric around the carrier so that the quadrature processes 
are independent of each other. This assumption is valid for 
various fading spectra. For example, see [11] and its refer- 
ences. However, in some fading environments such as the land 
mobile channel with Jakes model [12], the Doppler spectrum 
becomes asymmetric when the multipath signals are absorbed 
by obstacles or the propagation environment is characterised 
by directional non-isotropic scattering [13] — [15]. Thus, it is 
of great practical importance to take account of the effect of 
the asymmetric fading spectrum on the receiver structure and 



the performance analysis of differentially detected DPSK over 
i.n.i.d. channels, the topic of this paper. 

The paper is orgainzed as follows. In Section II, the signal 
model is introduced and different optimum diversity receivers 
are derived for different Rayleigh fading scenarios (see eqs. 
(17)-(20) below). In Section HL we use 8-DPSK as an 
example to study the BEP performance. Here, the average BEP 
of the optimum diversity receiver is obtained by calculating the 
BEP for each of the three individual bits. The results are given 
in exact, explicit, closed-form expressions which show clearly 
the behavior of the performance as a function of signal-to- 
noise ratio (SNR), fading correlation coefficient, and diversity 
order. Section IV presents numerical examples. Throughout 
this paper, overhead ~ denotes a complex quantity, superscript 
* will denote its conjugate, E is the ensemble average operator, 
5 represents the Kronecker delta, and [ ] T denotes transposition 
of the vector and matrix. 

II. SIGNAL MODEL AND RECEIVER STRUCTURE 

With space diversity reception over L frequency nonse- 
lective, i.n.i.d., Rayleigh fading branches with AWGN, the 
received signal over the ith branch, i = 1,2, ••• ,L, during 
the fcth symbol interval kT < t < (fc + 1)T is given, after 
matched filtering and sampling at time t = (k + 1)T, by the 
statistic r.i(k), where 



f i {k)=E 1 J 2 e^c l {k)+n l {k). 



(1) 



Here, E s is the energy per symbol, and for DPSK, 0(fc) is 
the data-modulated phase with Gray encoding of bits onto 
the phase transition A0(fc) = 4>(k) — (f>(k — 1). The fcth data 
symbol is conveyed in A0(fc). We assume here that all symbol 
points are equally likely. In (1), a rectangular data pulse shape 
g(t), where g(t) = 1/VT for < t < T and zero elsewhere, 
is assumed so that each matched filter has a rectangular low- 
pass-equivalent impulse response = g(T — t) for all i. 
Thus, the filtered noise ri;(fc) is given by 

^ T fti(t) 



hiik) 



kT 



T 



dt. 



(2) 



Here, is a set of i.n.i.d., lowpass, complex AWGN 

processes with E [fii(t)] = and E[ni(t)n*(t — t)] = iV^r) 
so that {ni(k)} k is a sequence of zero-mean, complex Gaus- 
sian variables with covariance function for each branch i 



E[n i (k)h* i (j)]=N i S kj 
The multiplicative distortion 6i(fc) in (1) is given by 

[ {k+l)T m 

Lt T 



Ci(k) 



dt. 



(3) 



(4) 



Here, {cj(i) = di{t) + J&i(t)}._ 1 is a set of i.n.i.d., lowpass, 
zero-mean, stationary, complex, Gaussian random processes. 
Each Ci(t) represents the complex gain due to frequency 
nonselective Rayleigh fading of the ith branch. For asymmetric 
spectrum in each i, the inphase fading process a.j(-) and the 
quadrature phase fading process fej(-) are generally correlated. 
At any time instant t, however, dj(i) and are always 

uncorrected. With reference to Fig. 1, it is shown in [16] that 
the covariance function of aj(-) and &,(•) can be obtained as 

E[ai(t)bi(t)] = (5a) 
E[ ai (t - T)oi(*)] = E[bi(t - T)bi(t)} = Ri(r) (5b) 
E[ ai (t)bi{t - r)] = -E[bi(t) ai {t - t)] = Qi(r) (5c) 

Note that if the spectrum of each Cj(i) is symmetric, the 
processes Oj(-) and fej(-) will be independent (i.e., we have 
Qi{r) = 0) with the same covariance function Ri(r). 

Letting c.;(fc) = ct;(fc) + jbi(k), it follows from (4) and (5) 
that both {di(k)}k and are sequences of zero-mean, 

real-valued, Gaussian random variables with 



E[oi(k)bi(k)] 
E[oj(fc — J)cij(fc)] 







E[6i(fc - i)6i(fc)] = Q(0 

(fc+l)T r{k+l-l)T 



II 



V) 



T 2 



kT J(k-l)T 

E[ ai (k)h(k - 0] = -E[6i(fe)oi(Ai - /)] = Di(l) 

kT J(k-l)T 

Thus, the covariance matrix can be obtained as 



J^2 



(6a) 
(6b) 

du dv 

(6c) 

du dv 



T, = E 



Oj(A; - I) 

bi(k) 
bi{k-l) 

Ci(0) 

Ci(l) Ci(0) 



a* (A;) ai(fc - I) bi(k) h(k - I) 



(7) 



A(0 

-A(0 o 
o -AO) Ci(o) Ci(i) 

AG) CiG) Ci(0) 

For each i, c;(fc) and nj(fc) are mutually independent. For 
* 7^ 3i {ci(k) , fii(k)} are independent of {cj(k),hj(k)}. 
The diversity branches are nonidentical since the covariance 
functions Ri(r), Qi(r) and NiS(r) depend on the branch 
index i. For convenience of later application, the following 
parameters are defined. The fading correlation coefficient at 
the matched filter output over a symbol interval of T for the 
ith diversity branch is defined as 

E[ci(k)ct(k-1)] Ci(l)-jDi(l) 



Pi 



EpCiWlVEOcitfc 



1] 



cm 



From (8), we note that pi is a complex quantity. It is a measure 
of the fluctuation rate of the channel fading process. The mean 
received SNR per symbol over the ith branch is defined as 



E[\El /2 e^c t (k)\ 2 ] _ 2E s Cm 



(9) 



Ni Ni 
We consider 2-, 4- and 8-DPSK with Gray encoding of bits 
onto A(j)(k) as shown in [4, Fig.l] for 4- and 8-DPSK, the 
mean received SNR per bit 7 £ is given by 7 f = 7i for 2-DPSK, 
7 f = 7t /2 for 4-DPSK, and 7l & = 7 ,/3 for 8-DPSK. 



Using the MAP criterion, the aim of the receiver is to 
determine from the received signals {f,(fe), fi(k — l)}f = i 
which one of the possible values 2-KmfM, m = 0, 1, • ■ ■ , M— 
1, of the phase difference A<fi(k) has maximum probability 
of occurrence. Following [9], it can be shown that MAP 
detection is equivalent to maximum log-likelihood detection. 
Specifically, based on {fi(k), fj(fc — l)}f =1 , we decide that 
A<fi(k) = 2nn/M whenever the log-likelihood function 

L 



iog* m = lo g \p 



h{k) h{k-l),A<j>{k) 



2irm 



M 



(10) 



is maximized for m = n. 

To proceed with evaluating (10), we need to verify that 
Ci(k) = a,i(k)+jbi(k) andci(fc— 1) = a,(fc— l)+jbi(k— 1) are 
jointly complex Gaussian. By being jointly complex Gaussian, 
it means that if x = xr + jxi and y = yR + jyi are two 
column complex random vector, then [xr t yR T xi T yi T ] T 
has a real multivatiate Gaussian probability density function 
(PDF), and furthermore, if u = [xr t yR T ] T and v = 
[xi T yi T ] T , then the real covariance matrix of [u T v T ] T 
has a special form given in [18, Theorem 15.1] that satis- 
fies Goodman's theorem [19]. After careful examination, it 
follows from (7) that Cj(fc) and 5j(fc — 1) are indeed jointly 
complex Gaussian^. Thus, conditioned on Ci(k — 1), (k(k) is 
conditionally complex Gaussian with mean [18] 



E[ci(k)\ci(k-l)]=piCi(k-l) 



(ID 



and variance 
E 



{|ci(*)~ 
= 2Ci(0) 



E[ci(*)|ci(A - l)]| 2 |ci(A; - 1)} 
Cf(l)+DKl) 



(12) 



Moreover, conditioned on the vector [<2j(fc — 1) bi(k — 1)] T , 
the vector [dj(fc) bi(k)] T is conditionally Gaussian with 
covariance matrix given by 



a- 



Ci(0) 



g, (i)+-p, (i) 
cm 

o 



cm - 



o 

g, 2 (i)+-p?(i) 
cm 



(13) 

which is a diagonal matrix. This shows that Ke[ci(k)\ci(k— 1)] 
and lm[ci(k)\ci(k — 1)] are independent. 

Applying (11), (12) and (13) to (10), we obtain 



^log*r 



Re 



E 



= C+ (14) 

2E S [d(l) + jA(l)] e-J 2 ^ r,(fc)f*(fc - 1) 
[2E S CM + Ntf - 4A 2 [C?(1) + A 2 (l)] 



(8) or, equivalently 



(15) 



Re 



y- 



fi(fc) f*(fc - l)e 3 ™ 



where £ represents the constant term whic h does not af fect 

/ C 2 (l)-\-D 2 (1) 

the decision. In (15), the quantities \pi\ = J ' an & 



Zp.i = — tan 



-l 



cI(TT 



represent the magnitude and phase of 



the correlation coefficient pi given in (8), respectively. 

'We also call them the proper complex Gaussian random variables [17]. 



Defining the real-valued weighting factors 

J \pi\li 

Ni 



(16) 



(l + 7 i ) 2 -(|p,l7 i ) 
it follows from (15) that the optimum combining differential 
receiver will now compute, for the fcth symbol, the decision 
statistics {A m (fc)}m=o' and declares that Acj>{k) = 2™ if 
A„(fc) = max m {A m (fc)}, where 

L 

A m (fc) = Rc w i f i( k ) f i ( k ~ !) e 



(17) 



If the spectrum of the channel complex gain is symmetric, 
Pi is a real-valued quantity. Then, the optimum combining 
differential receiver (17) will become [9] 

L 

A' m (k) = Re e~ j ^ ^w, h{k)f*{k - 1) (18) 

i=l 

If the diversity branches are i.i.d., but the fading gains have 
asymmetric spectrum, the optimum receiver will become 



Al(fc) = Re 



(19) 



where p — pi for i — 1,2,- • • , L. For i.i.d. branches with 
fading gains having symmetric spectrum, the optimum receiver 
is the well-known product detector, given by [4] 

L 



AZ(k) = Rc 



(20) 



Comparing (20) with (17), we see that in the case of i.n.i.d. 
channels with asymmetric power spectrum, the receiver first 
rotates the product phasor fi(k)f*(k — 1) between the two 
received signal samples at each diversity branch by the angle 
— /.pi, then scales each resulting phasor by the weight Wi, and 
finally sums all L rotated and scaled phasors to form a decision 
variable. Clearly, in order to form the optimum detector (17), 
besides the received signal samples fi(k) and fi(k — 1), 
the receiver requires the a priori knowledge of the channel 
statistics, including the power spectral densities of AWGN 
Ni, both the magnitude and phase of the fading correlation 
coefficient pi, and the mean received SNR 7.;. These quantities 
can be pre-computed according to our knowledge of the 
channel statistics at the receiver. 

III. PERFORMANCE ANALYSIS 

In this section, we will derive exact, explicit and closed- 
form BEP expressions for differentially detected DPSK for 
the optimum receiver (17). Due to space limitation, we only 
consider 8-DPSK in this paper. The signal constellation, bit 
mapping and the decision region R m for 8-DPSK is shown in 
Fig. 2. In [4] and [9], the average BEP is computed using 
the binary reflected Gray code (BRGC) approach through 
Hamming weight spectrum [20]. It is shown in [21] that the 
BRGC approach with Hamming weight is less accurate for 
M > 16. In this paper, we adopt a new approach, namely, 
the average BEP is obtained by calculating the BEP for each 
of the three individual bits in 8-DPSK. This approach has 
the advantage of showing explicitly the BEP performance 
differently for the three different transmitted information bits. 
Therefore, using the bit which has lower BEP to convey more 
important information can improve communication reliability. 



From Fig. 2, we see that each signal point is represented 
by a 3-bit symbol (ji, j 2 , j'3). We use Pj 1 , Pj 2 and Pj 3 to 
denote the corresponding individual BEP. Since the three bits 
are equally likely, the average BEP is given by 



P = -(P h 
3 V Jl 



Pi. 



(21) 



We begin with computing Pj 1 . Without loss of generality, it 
is assumed that j\ = 0. The case where j\ = 1 gives an 
identical result. From Fig. 2, we see that the bit ji = is 
associated with the symbols 000 (A</>(fc) = 0), 001 (A</>(fc) = 
tt/4), 011 (A(j)(k) = tt/2), and 010 (A<j)(k) = 3tt/4). Thus, 
conditioning on j\ = 0, the BEP Pj 1 will be given by 

1 \P h (e| Acj>{k) = 0) + P n (e| A<f>(k) = tt/4) (22) 



Pjl ~ 4 



+P n (e| A<f>(k) = tt/2) + P h (e\Ac/>(k) = 3tt/4) 



Here, P n (e\A(/>(k) = mix /A), m = 0, 1, 2, 3, is the probabil- 
ity that conditioning on A(f>(k) = mix /4, the decision j\ = 1 
is made. With reference to Fig. 2, this is equivalent to the 
probability that conditioning on A^(fc) = mn/A, the phasor 
^^ = iWifi{k)f*{k — X)e~i^ i lies outside the half-plane 
region R0+R1+R2+R3 (i.e., in the region R4+R 5 +Rq+Ri). 
The BEP P h (e\A(f)(k) = nnr/4) is thus obtained as 



r/4) = P\ Re 



(23) 




Wifi{k)f*{k - l)e 



< 



A0(fc) = — 



To evaluate (23), first, it follows from (11) and (12) that con- 
ditioning on A(f>(k) = mir/4: and on f.;(fc — 1) eJ z P { eJ 37T ^ 8 = 
a,, for i = 1,2, the quantity fi(k) is condition- 

i e -j37r/8 e J 



ally Gaussian with mean on y-p^- e 



Pili p -jZpi -j3tt/8 jmir/4: 



% e -l^/8 elm ^/4^ where 



Pi 



Pi\e Pi has been used, 



(l+7») 2 -(|p,|7>) 2 
l+7i 



Ni. Then, in (23) the quantity 



1+7, 
and variance 

Refe -5 ^ Wifi(k)f*(k — l)e~^^*)] is conditionally 
Gaussian with mean cos (m7r/4 — 37r/8) Yn=i w i T^fr" l^^ 2 

N t \& t \ 2 . Finally, 



ls^ L ,,,2 (1+7.) -(Ipih.r at 1^. 12 

2 Z^i=l w i l+ 7l 



and variance 

following the derivation procedure detailed in [9], the BEP 
in (23) can be obtained as 



P, i (e|A^(fc) = ^) = 



(24) 



L „ 



i=l 



1 



cos z 



< 4 



8 , 



cos- 



m7r 

> 4 



37T 

8 



) + 1/Ai 



where the quantity Gi is given by 

G l = l —, and A.; = 



3 = l,j^i 



A, — A, 



£ + li) 2 -{\~Pi\lif 



(25) 



Putting (24) into (22) leads to the BEP P jx . An interesting 
observation from (24) is that the BEP does not depend on the 
phase, Zpi, of the fading correlation coefficient pi. Intuitively, 
this is because the optimum receiver (17) can provide "phase 
compensation" for each diversity branch before combining 
using the channel statistic knowledge e~ 3 As such, we 
expect that the receivers (18) and (20) are suboptimum over 
the channel with asymmetric fading spectrum. 



Next, we compute Pj2 in (21). The procedure for obtaining 
the conditional BEP for j'2 = is parallel to that followed in 
the case for j\ = 0. From Fig. 2, the bit J2 = is associated 
with the symbols 001 (A<p(k) = tt/4), 000 (A<f>(k) = 0), 
100 (A0(fc) = 7n/A), and 101 (A<j>(k) = 3tt/2). Hence, 
conditioning on j% = 0, the BEP Pj 2 is given by 

^ = \ [Ph (e| A0(fc) = tt/4) + P J2 (e|A0(fc) = 0) (26) 

+P, 2 (e|A0(fc) = 7tt/4) + P, 2 (e|A0(fc) = 3tt/2) 

where P J2 (e| A<f>(k) = nn/A),n = 0, 1, 6, 7, is the conditional 
probability that the phasor Ylf=i w i^i(k)r* (k — l)e~ jZ ^ i lies 
in the half-plane region R2 + R3 + R4 + P5, i.e., 



P j2 (e\A0(k) 



r/4) = P< Rc 



(27) 



Vi=i 

which has solution 



l)e 



< 



A0(fc) = — 



P i2 (e|A0(fc) = m/4)= 



(28) 



T 

i=l 



1 



COS" 



' n7r 7r \ 
k 4 ' S) 



C:os- 



' n7r 
^ 4 



Putting (28) into (26) leads to the BEP P j2 . 

Finally, we compute Pj^ in (21). From Fig. 2, the bit J3 = 
is associated with the symbols 100 (A<f>(k) = 7n/4), 000 
(Acj){k) = 0), 010 (A0(fc) = 3tt/4), and 110 (A<f>(k) = n). 
Thus, conditioning on j'3 = 0, the BEP Pj 3 is given by 

Ph = \ [ p n (e| A0(fc) = 77T/4) + P J3 (e|A#fc) = 0) (29) 

+P is (e|A#fc) = 3n/4) + P h (e\A<f>(k) = tt) 

where P, 3 (e|A0(fc) = Z7r/4) , Z = 0,3,4,7, is the conditional 
probability that the phasor X^ =1 Wifi(k)f*(k — l)e~^^ lies 
in the region Pi + R2 + R5 + Re- This is equivalent to the 
conditional probability that after rotating by — ir/8, the product 
of the inphase and quadrature-phase components of the phasor 
Ylt=i w i^i(k)r* (k — l)e~-' z ' 5i is greater than zero, i.e., 

In/ A) = (30) 



P, 3 (e|A^(fc) 



P{ Re 



y^it;^(fc)f*(fc - l)e 

V.i=l 



Im 



l)e 



> 



A^(fc) = ^ 



x ^2wifi(k)f*(k 

\i=l 

From the argument for deriving (24), we note that condi- 
tioning on A<j>(k) = /tt/4 and on f t (k - 1) e j/ « e^/ 8 = 
/3i, for i = 1,2, the inphase component in (30), 

Re[e~ J f (X^i Wifi{k)f*{k — l)e _jZ ^ i )] is conditionally 
Gaussian with mean cos (In/ A — 7r/8) |^| 2 
w 2 (i+7,) 2 -(Ip,|7,) 2 Ar . 1/5.12 



and variance | 2»=i 



l+7i 

Similarly, 



the component Im[e _: '?'(^^ 1 Wifi(k)f*(k — l)e _J w )] in 
(30) is also a conditionally Gaussian random variable, 



with mean sin (In/ A 



2 



2 ^i=l W » 1+7, 

lows from (13) and the properties of the complex Gaussian ran- 
dom variables [18] that the conditional inphase and quadrature- 
phase components Re[e _:, t (^2f =1 uiifi(k)r* (k — l)e~ jZ ^)] 



2 (i+7.) -(ip,i7.r 



1+7, 



|/3i| 2 and vari- 
Ni\pi\ 2 . Moreover, it fol- 



and Im[e _:, t (5^ i=1 Wifi(k)f*(k — l)e _ ' jZ ^ i )] in (30) are also 
independent. Therefore, conditioning on A<j)(k) = In /A and 
on fi(k — 1) e- jZ ^ i e^/ 6 = Pi, and denoting the inphase and 
quadrature-phase components as 



X 
Y 



n(cos (In I A - 7r/8) u, rj- 
JvYsin (/7r/4 - n/8)u, rf 



(31) 



where u and 77 are given, respectively, by 

L 

u= 

L ■ \2 (\7. U, s |2 



\ - PiTh |„ 



E 



2 (l+Ti) 2 



1 



2^ l + 7l 



id: 



(\Pi\7iY 



(32) 



the conditional BEP P js (e A^(fc) = -f > A) is g iven b Y 

P i3 (e|A^(fc) = ^A)=P^y> O|A0(fc) = 

(33) 

This is probability that the product of two independent real- 
valued Gaussian random variables with non-zero, nonidentical 
means and identical variances is greater than zero. This is a 
special case of the results given in [2, Appendix B] concerning 
the probability that a general quadratic form in complex-valued 
Gaussian random variables is less than zero. Using [2, (B-21) 
of Appendix B], (33) can be evaluated as 



P i3 (e|A0(fc) = ^,A) =1 



(34) 



Qi (V.g[l-sin (Ztt/2 - tt/4)], ^ g[l + sin (In/ 2 - tt/4)]) 



1 



I Q [g\ cos (Ztt/2 - ir/A)\] exp(-g) 



where, Qi(o, 6) is first-order Marcum's Q-function and Jfc(x) 
is the fcth-order modified Bessel function of the first kind. In 
(34), the quantity g = Yd=i K \h? has PDF given by [9] 

L 



-<^(l + 7l ) 



exp 



9 



w' l N l (l+ ll ) 



(35) 



where w' ; = -=V 



\T7t^t^tI^ tt= 1 — rar- Averaging the condi- 

Ni (l+7i)[(l+7i) 2 -(|Pi|7i)' ! ] & & 



tional probability (34) over g using the PDF (35) gives the 
BEP P js (e\A<j)(k) = In /A) in (30), i.e., 

In f°° In ~ 

P h (e\A<t>(k) = —) = J P h (e\A<f>(k) = -,Pi)p(g)dg 

(36) 

Substituting (34) and (35) into (36), we obtain, after manipu- 
lation and simplification, 

L 



p J3 ( e \Am = l i) = Y< G 



A, 



1 



1 



|cos(i 



: )I(V2-1) 



where A; is given by 



A\ - cos2 ( l f - f ) 
1/2 

A 2_ cos 2(i f _ f)J 



(37) 



1 + 7. 



Putting (37) into (29) leads to the BEP P j3 . Substituting (22), 
(26) and (29) in (21), we obtain the average BEP P. 



IV. NUMERICAL EXAMPLE 

Fig. 3 plots the BEP performance for the three individual 
bits in (22), (26) and (29) and the average BEP in (21) of 
8-DPSK against the total average received SNR per bit. The 
order of diversity is set to L = 2. The abscissa represents the 
total mean SNR per bit which is given by 7& = Y2i=i l\ = 
i Y^i=i li- The average received bit energy distribution among 
the two branches is set to : ~j\ = 30% : 70%. It is 
assumed that the fading correlation coefficient (the normalized 
covariance function) model follows [14, eq.(10)], given by 



Efc(*)c?(t-T)] 
E[|Si(*)| 2 ] 



lof^y/^-A^f^+jAnKfar) 



(38) 

where fa is the Doppler frequency, and k is a parameter that 
controls the width of the angle of arrival of scatter components 
[14, eq.(l)]. Note that if k = 0, (38) results in the correlation 
coefficient for the Jakes two-dimensional isotropic scattering 
model, i.e., E[ci(t)%(t - r)]/E[|ci(t)| 2 ] = I (j2irf m T) = 
Jo(27r/ m r), where Jo(-) is the zeroth-order Bessel function. 
We assume that the normalized Doppler spread faT = 0.03 
and 0.05 for diversity branches 1 and 2, respectively, and the 
parameter k is set to 3. Thus, we have p\ = 0.9871 +j0. 1519 
and p 2 = 0.9642 + jO.2511. It is seen from Fig. 3 that the 
third bit j'3 has the lowest BEP, whereas, the BEP Pj 1 for the 
first bit ji is equal to the BEP Pj 2 for the second bit j 2 . 
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Fig. 1. Illustration of complex channel fading process. 
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Fig. 2. 8-DPSK constellation and decision region. 
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Fig. 3. BEP comparison of the three individual bits and the average of all 
bits for 8-DPSK. 



